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This  report  presents  results  of  work  performed  by  the  Lockheed- 
Huntsville  Research  & Engineering  Center  under  Contract  DAAH03-75-A-0040, 
Delivery  Order  CC-10,  for  the  U.S.  Army  Missile  Command,  Redstone  Arsenal, 
Alabama.  From  the  administrative  point  of  view,  this  document  represents  a 
final  report.  However,  the  work  itself  has  raised  a sufficient  number  of  as  yet 
unanswered  questions  such  that,  from  the  technical  point  of  view,  this  should 
only  be  considered  as  an  interim  report. 

This  work  was  monitored  by  Mr.  Charles  M.  Cason  of  the  Laser  Science 
Division,  Army  High  Energy  Laser  Directorate.  The  period  of  performance 
covered  by  this  report  was  from  11  December  1975  through  30  August  1976. 
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Section  1 

INTRODUCTION  AND  SUMMARY 

Closed  cycle  operation  of  CO^  electric  discharge  lasers  has  been 
shown  to  be  substantially  influenced  by  the  buildup  of  molecular  species  with 
high  electron  attachment  cross  sections,  negative  ions  which  increase  the 
plasma  heating,  and  species  which  efficiently  deactivate  the  upper  laser  level. 
The  plasma  chemistry  important  in  the  formation  of  these  impurities  has  been 
studied  in  greatest  detail  for  He/N^/CO^  mixtures  loaded  by  a dc  discharge. 
Current  Army  emphasis  is  on  elimination  of  helium  containing  mixtures  (for 
logistical  reasons)  and  the  development  of  repetitively  pulsed  laser  technology. 
Thus,  particularly  for  closed  cycle  operation,  there  exists  interest  in  deter- 
mining dominant  processes  detrimental  to  sustained  repetitive  high  energy 
density  pulsed  loading  of  generic  H^/N^/CO^  laser  mixtures  which  might  also 
contain  CO,  H^O  or  other  species  characteristic  of  hydrocarbon  fuel  com- 
bustion products. 

The  purpose  of  the  effort  described  here  was  to  identify  chemical 
reactions  leading  to  impurity  production  in  several  laser  gas  mixtures, 
to  determine  the  relative  importance  of  the  impurity  reaction  products  in  de- 
grading assessed  laser  performance,  and  to  furnish  predictions  of  all  im- 
portant molecular  species  production  and  decay  rates  during  and  following 
the  electrical  pulse.  In  doing  this,  maximum  use  was  to  be  made  of  existing 
kinetic  data  and  of  the  existing  continuous  wave  electric  laser  code  previously 
developed  by  the  authors  (Ref.  1). 

Furthermore,  the  steady  state  buildup  of  impurity  concentrations  under 
repetitively  pulsed  operation  was  to  be  examined,  and  key  experimental  efforts 
to  verify  predicted  results  were  to  be  recommended. 
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The  essential  features  of  the  modeling  approach  used  in  this  effort  are 
described  in  Section  2.  The  discussion  includes  the  basic  equations  as  well 


I 


as  the  details  of  the  plasma  kinetics  model  employed.  Results  obtained  for 
several  laser  gas  mixtures  are  described  in  Section  3.  Discussion  of  the 
standard  CO^/N^/He  laser  gas  mixture  includes  a comparison  of  performance 
characteristics  as  predicted  by  the  model  developed  under  this  contract  with 
experimental  data.  The  good  agreement  achieved  lends  credibility  to  the  theo- 
retical model  and  its  application  to  the  analysis  of  other  laser  gas  mixtures. 

Section  4 presents  the  conclusions  derived  from  the  work  accomplished. 
Also  presented  are  recommendations  for  key  experimental  efforts  that  should  be 
undertaken  to  verify  predicted  results,  as  well  as  for  future  theoretical  efforts 
that  should  be  performed  to  develop  a more  user -oriented,  efficient  computer 
model  specifically  suited  to  investigate  pulsed  electric  lasers. 

Briefly  summarizing,  a theoretical  model  is  presented  which  can  be  used 
to  analyze  plasma  chemistry  effects  on  the  performance  of  pulsed,  electron  beam 
sustained  electric  discharge  lasers.  The  one-dimensional  model  contains  de- 
scriptions of  the  configuration  components  of  primary  importance  which  are  the 
electron  beam,  the  discharge,  the  plasma  and  its  composition,  laser  radiation, 
and  their  mutual  interaction.  All  important  processes  are  modeled  in  terms 
of  experimentally  controllable  parameters  which  include  the  post-foil  electron 
beam  current  density,  the  primary  electron  energy,  the  applied  electric  field 
potential,  discharge  electrode  spacing,  initial  laser  gas  composition,  tempera- 
ture and  pressure,  both  electron  beam  and  electric  discharge  pulse  widths  and 
relative  sequencing,  as  well  as  cavity  geometric  data  such  as  mirror  reflec- 
tivities, transmissivities  and  separation. 

Calculated  single  pulse  performance  for  a CO^  laser  oscillator  at  5 atm 
is  found  to  be  in  good  agreement  with  experiment.  Analysis  of  several  potential 
laser  gas  mixtures  containing  hydrogen  in  place  of  helium,  under  conditions 
representative  for  the  small  scale  closed  cy  . > circulator  planned  by  the  Army, 
indicated  that  closed  cycle  operation  and  repetitive  pulsing  leads  to  a strong 
degradation  of  laser  performance  due  to  formation  of  electrophilic  species. 
Important  contaminants  are  identified,  and  their  selective  removal  is  shown 
to  prevent  degradation  of  laser  performance. 
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2.1  THEORETICAL  MODEL 

This  section  presents  a brief  description  of  the  essential  features  which 
a theoretical  model  must  possess  to  adequately  describe  plasma  chemistry 
processes  in  pulsed  electric  discharge  lasers.  This  model  incorporates  de- 
scriptions of  the  electron  beam,  the  discharge,  the  plasma  composition  and 
laser  radiation,  and  their  mutual  interaction.  The  effect  of  both  single  and 
repetitive  pulsing  can  be  modeled.  To  simulate  selective  removal  of  detri- 
mental species  (scrubbing)  the  plasma  composition  can  be  modified  between 
successive  pulses  when  recirculation  and  repetitive  pulsing  are  to  be  con- 
sidered. Modeling  as  many  of  the  basic  features  of  a device  as  possible  serves 
to  minimize  simplifying  assumptions  which  always  have  to  be  made.  An  impor- 
tant simplification  made  in  this  work  is  the  assumption  of  spatially  uniform 
conditions  in  the  cavity  resulting  in  a one-dimensional  model.  On  the  other 
hand,  the  full  coupling  of  the  energy  equation  into  the  system  of  species  rate 
equations  eliminates  the  often  made  assumption  (see  Section  4)  that  variations 
in  the  gas  temperature,  which  is  of  prime  importance  to  laser  operation,  can 
be  ignored. 

Considering  uniformly  excited  laser  cavity  volumes  of  the  order  of  a few 
liters,  and  electric  discharge  pulses  of  approximately  microsecond  duration  it 
can  easily  be  shown  that  a constant  density,  time -dependent  treatment  of  the 
plasma  chemistry  should  adequately  describe  the  important  processes  occur- 
ring in  a pulsed  discharge.  For  repetitively  pulsed  discharges  and  closed  loop 
operation  of  the  laser,  the  flow  velocity  around  the  cycle  is  assumed  to  be  suf- 
ficiently small  that  pressure  changes  due  to  varying  flow  velocity  can  be  neg- 
lected. Plasma  composition  changes  around  the  recirculation  loop  can  then  be 
treated  in  the  same  time-dependent  manner  as  in  the  cavity,  provided  that  post- 
discharge expansion  of  the  plasma  as  well  as  cooling  by  heat  exchangers  is 
appropriately  accounted  for. 
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The  laser  gas  mixture  generally  consists  of  neutral  atoms  and  molecules, 
charged  atoms  and  molecules,  vibrationally  excited  molecules,  and  electrons. 
Therefore  a multitude  of  species  must  be  tracked  during  the  calculations.  It  is 
assumed  that  different  quantum  states  of  various  particles  can  be  treated  as 
separate  species.  For  simplicity,  it  is  also  assumed  that  a single  vibrational- 
rotational  radiative  transition  sufficiently  accounts  for  laser  output.  While  not 
generally  true,  this  is  a reasonable  assumption  if  lasing  of  the  CO^  molecule  is 
consider ed. 

Under  the  above  stated  assumptions  the  basic  governing  equations  can 
be  written  as  follows: 


Continuity: 


Species: 


(1) 


(2) 


Here  p denotes  the  density,  Fi  the  mol/mass  ratio  of  the  ith  species,  and  wi 
its  overall  net  rate  of  production  due  to  any  reaction  mechanism  except  for 
radiation.  For  an  arbitrary  reaction  of  the  form 


the  net  rate  of  production  for  any  participating  species  for  which  i/!'  - i/!  / 0 
can  be  written  as 


( N I'l  N ,,) 

K - •'P  ktn,  »>FP  -kb  n mpv 

< i=l  i=l  I 


where  kj  and  k^  denote  the  forward  and  backward  rate  constant,  respectively. 
Laser  radiation,  its  intensity  denoted  by  I,  depletes  the  upper  laser  level  and 
replenishes  the  lower  laser  level,  as  indicated  by  the  Kronecker  symbols  in 
Eq.  (2).  Gain  and  photon  energy  are  denoted  by  a and  e,  respectively. 
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Energy: 


£(pe»=  jnE-„aI 


which  simply  states  that  any  change  of  internal  energy  of  the  laser  gas  is  due 
to  discharge  energy  pumped  into  the  plasma,  or  due  to  radiation  leaving  the 
gas  in  the  cavity.  Here  jj^  denotes  the  drift  field  current  density  and  E the 
electric  field  strength.  The  quantity  K represents  a correction  factor  and  is 


defined  as 


K s 1 - exp(-qL) 


where  a is  the  cavity  threshold  gain.  Assuming  a Fabry-Perot  cavity,  the 
threshold  gain  can  be  expressed  in  terms  of  mirror  reflectivities  and  optical 
path  length  through  the  gain  medium  as 


a = - 


fn(rir2) 


Note  that  for  aL  « 1,  which  is  most  often  the  case,  k reduces  to  unity.  Lasing 
commences  when  the  medium  gain  exceeds  the  cavity  gain,  and  ceases  when  the 
pumping  reactions  fail  to  maintain  the  medium  gain  at  the  cavity  gain  level. 

The  medium  gain  can  be  expressed  as 


a = p € ( S 


F 

21*2 


S12Fi)V 


where  F^  and  Fj  are  the  mol/mass  ratios  of  the  upper  and  the  lower  laser 
level,  respectively.  All  spectroscopic  data  enter  the  equation  through 
and  Sj2  which  are  functions  of  temperature  only.  V denotes  the  Voigt  function 
which  basically  represents  a generalized  line  shape  function, 


V(rj)  = exp(rT)  [l  - erf(n)] 
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with  rj  being  the  ratio  of  the  Lorentz  line  width  to  the  Doppler  line  width,  i.e., 


r)  = (in  2) 


1/2  AvL 


Avt 


(10) 


The  set  of  basic  equations  is  completed  by  the  equation  of  state, 

NS 

P = P T J F.  (11) 


i = l 


which  determines  the  pressure  as  a function  of  density,  temperature  and  gas 
composition. 


Energizing  the  laser  gas  in  an  electric  discharge  requires  the  presence 
of  free  electrons  which  are,  at  least  initially,  here  assumed  to  be  generated  by 
electron -beam  ionization.  Assuming  that  both  the  effective  electron-beam 
current  density,  jgg,  and  the  effective  electron  beam  energy,  Egg,  are  speci- 
fied, the  number  density  of  primary  electrons  in  the  laser  gas  can  be  expressed 
as 


n 


JEB 


EB 


q v 

Me  vEB 


(12) 


where  is  the  electron  charge,  and  Vgg  the  beam  electron  velocity  expressed 
as 

1/2 


EB 


= c 


/2Eeb\ 

\eref/ 


(13) 


Here  c denotes  the  speed  of  light,  and  E^gg  the  electron  rest  mass  energy 


equivalence  (mg  c = 511  keV).  Knowing  the  primary  electron  density,  the 


rate  of  ionization  and  production  of  secondary  electrons  due  to  the  electron 
beam  only  can  be  determined  according  to  the  reaction 


Ai  + nEB 


Li  .+  , 

Ai  + ne  + nEB 


(14) 


Hi 

i 


where  the  primary  ionization  rate  constant,  can  be  evaluated  as  a func- 

tion of  the  ion  mass,  the  beam  ionizing  power  and  the  beam  electron  velocity 
(Ref.  1). 


The  drift  field  current  density,  and  therefore  the  energy  which  can  be 
pumped  into  the  gas,  is  proportional  to  the  number  of  free  electrons  and  their 
drift  velocity  which  in  turn  is  dependent  on  the  electrical  field  strength,  such 
that,  to  good  approximation, 


J 


D = C ne  qe  (E/N) 


3/4 


(15) 


The  volumetric  rate  of  energy  addition  in  the  discharge  can  then  be  expressed 
as 


jDE  = CneqeN(E/N) 


(16) 


where  the  constant  C,  strictly  speaking,  is  mixture-dependent,  but  in  practice 
was  found  to  be  essentially  independent  of  the  gas  mixture  over  the  range  of 
compositions  considered.  Equation  (16)  is  not  only  used  to  evaluate  the  energy 
added  in  the  discharge  but  also  to  determine  the  secondary  electron  energy 
which  is  essential  in  evaluating  electron  impact  excitation  rates.  Postulating 
a basic  electron  impact  vibrational  excitation  mechanism  as 

k'  , 

(17) 


f,i 


A^(v)  + e A^v'  > v)  + e 


the  power  consumed  by  this  mechanism  can  be  expressed  as 

**  = £ “e  ”i  kf,  i iHi  <‘S> 

1 

where  the  summation  indicates  that  more  than  one  mechanism  may  be  involved. 

Applying  conservation  of  energy  by  equating  the  right  hand  sides  of 
Eqs.  (16)  and  (18)  then  leads  to  the  important  relation 


7 


(19) 


CqeN(E/N)7/4  = k£  . AHj 


Since  the  impact  excitation  rate  constants,  k,  are  functions  of  the  electron 
energy  which  itself  is  a function  of  E/N,  Eq.(19)  constitutes  an  implicit  relation 
for  the  electron  energy,  and  simultaneously  states  its  functional  dependence  on 
E/N  in  terms  of  input  functions. 


Equations  (1)  through  (19),  in  conjunction  with  a detailed  set  of  reactions 
and  rate  coefficients,  to  be  discussed  in  Section  2.2,  provide  the  basic  frame- 
work to  evaluate  plasma  conditions  in  the  discharge  as  a function  of  time,  and 
if  desired,  the  radiative  output  of  the  laser  cavity. 


In  theoretically  assessing  the  effects  of  recirculation  on  electric  dis- 
charge laser  performance  it  is  most  important  to  achieve  a reasonable  simu- 
lation of  the  density  and  temperature-time  history  around  the  recirculation 
path,  because  the  reaction  rates  are  essentially  functions  of  these  two  param- 
eters. The  procedure  followed  in  the  present  investigation  was  to  pulse  the  gas 
in  the  laser  cavity  for  the  duration  of  microseconds  at  constant  density.  From 
this  state  of  elevated  pressure  and  temperature,  the  plasma  is  allowed  to  re- 
turn to  its  initial  pressure  via  an  adiabatic  expansion.  The  time  interval 
appropriate  for  this  expansion  process  can  be  estimated  by  considering  the 
time  which  a disturbance  traveling  at  the  speed  of  sound  would  take  to  clear 
the  laser  cavity  or  discharge  volume.  Typically,  this  should  occur  in  a time 
interval  of  the  order  of  100  microseconds.  The  theoretical  cycle  is  completed 
by  considering  an  isobaric  cooling  process,  simulating  the  flow  of  gas  through 
a heat  exchanger,  which  returns  the  gas  to  the  initial  temperature  (and  pres- 
sure). Typical  time  intervals  for  this  last  leg  of  the  cycle  would  be  on  the 
order  of  milliseconds. 

2.2  PLASMA  KINETICS 

In  Refs.  1 and  2 we  discussed  in  some  detail  the  input  kinetics  developed 
(under  Contracts  DAAH03-75-C-0077  and  DAAH03 -75-A-0040,  Delivery  Order 
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CC-ll)  to  compute  the  gasdynamic,  discharge  and  gain  characteristics  of  a 
CW  CO^  electron  beam-sustained  discharge  with  air  contaminant.  Operation 
of  the  laser  with  bottled  CO^-CO-N^-He  gas  mixtures  of  standard  purity  at 
a pressure  on  the  order  of  0.1  atm  was  implicitly  assumed.  Here  we  have 
extended  the  kinetics  to  consider  repetitively  pulsed  operation  at  pressures 
up  to  10  atm  with  cryogenically  dried  laser  mixtures  generated  by  combustion 
of  hydrocarbon  fuel  with  air.  In  this  section  we  present  the  updated  plasma 
chemistry  package  currently  in  use  (Tables  1,  2 and  3)  together  with  a very 
brief  discussion  of  the  changes  only. 

2.2.1  CO^-CO-N^-He  Mixtures  (Non-Hydrogenic) 

With  relatively  minor  modifications  the  plasma  chemistry  package  of 

Ref.  2 has  been  utilized  as  a logical  starting  point  for  similar  analysis  of 

closed  cycle  pulsed  operation  of  CO^  electric  discharge  lasers  utilizing  1 - 

2-3  CO^  - - He  mixtures  at  atmospheric  and  higher  pressures,  with  the 

very  encouraging  results  discussed  in  Section  3.  In  a sense,  the  apparently 

successful  extrapolation  of  these  basic  plasma  kinetics  to  a density  regime 

almost  two  orders  of  magnitude  higher  is  remarkable.  We  had  anticipated 

that  heavy  ion  clustering  (Ref.  3)  would  play  a more  dominant  role  and  that  as 

a minimum  revision  it  would  be  necessary  to  regard  the  rate  coefficients 

entered  for  the  parent  ion  of  a given  cluster  family  (e.g.,  CO*  • (CO_)  ) as  a 

u Ct  n 

global,  scalable  descriptor  for  the  rates  of  reaction  of  members  of  that  par- 
ticular cluster  ion  family.  Such  revisions  may  yet  be  required.  It  is,  how- 
ever, already  apparent  that  the  kinetics  code  provides  useful  predictions  using 
inputs  not  radically  revised  from  those  employed  in  modeling  plasmas  at  much 
lower  pressures.  To  rationalize  this  finding,  it  is  interesting  to  speculate  that 
the  heavily  clustered  ions  observed  in  the  recent  drift  tube  experiments  by 
McDaniel  et  al.,  (Ref.  3)  may  be  unstable  in  the  high  energy  density  field  of 
CO^  laser  discharges.  An  analogous  instability  of  complex  cluster  ions  in  the 
D region  of  the  ionosphere  apparently  occurs  in  the  presence  of  intense  solar 
radiation  during  solar  flares,  accounting  quantitatively  for  observed  enhanced 
free  electron  concentrations  (Ref.  4). 
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The  modifications  to  the  non -hydrogenic  kinetics  inputs  as  listed  in 
Ref.  2 are  as  follows: 

Ionized  Species  (Table  1) 

Reaction  B2;  In  addition  to  the  NO^  dissociative  detachment  reaction  B1 
theoretically  evaluated  by  Nighan  and  Wiegand  (Ref.  5),  the  associative  attach- 
ment process  B2  forming  NO^  has  now  been  incorporated.  Attachment  coeffi- 
cient B2  is  the  pressure-independent  experimental  value  measured  by  Mahan 
and  Walker  (Ref.  6)  in  the  presence  of  diluent. 

Reaction  Cl:  The  associative  detachment  coefficient  Cl  for  the  O"  + CO_ 

-0  4 -0  39 

reaction  has  been  trivially  rewritten  for  a T ’ rather  than  the  T ' tem- 
perature dependence  reported  by  McFarland  et  al  (Ref.  7). 

Reaction  D6:  Interchange  reaction  D6  between  CO"^  and  O atoms  has 

been  substituted  for  the  corresponding  process  involving  N atoms.  The  latter 

process  is  negligible  by  comparison,  as  is  readily  evident  from  the  lower  rate 

- 12 

coefficient  (viz.,  5.0  x 10  ) for  this  process  vis  a vis  that  for  D6,  each  of 

which  is  from  the  1973  survey  by  Ferguson  (Ref.  8). 

Reaction  D28:  Clustering  rate  coefficient  D28  has  been  changed  so  that  re- 
cently employed  by  Spjeldvik  and  Thorne  (Ref.  9)  in  their  D-region  ionic  model. 

Reaction  D29:  The  clustering  reaction  between  C^  and  CO^  with 
chaperone  species  has  been  incorporated  using  the  rate  coefficient  reported 
by  Fehsenfeld  and  Ferguson  (Ref.  10). 


Reactions  El,  E2  and  E3:  An  ion -temperature  dependence  has  been 

incorporated  into  the  very  important  ion-electron  recombination  coefficients, 

retaining  a u electron  energy  dependence.  In  the  absence  of  definitive 
- 1 G 

data,  a T ion  temperature  dependence  has  been  assumed  for  the  more 
complex  ions  (i.e.,  CO*  and  hydrated  ions),  this  being  the  limiting  depend- 
ence predicted  by  the  theoretical  analysis  of  Bates  (see  e.g.,  Ref.  11)  for 
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recombination  in  which  neutral  product  stabilization,  rather  than  radiationless 
transition,  is  rate -limiting.  Intuitively,  it  seems  that  the  multiplicity  of 
rovibronic  states  available  to  complex  ions  should  favor  radiationless  transi- 
tion by  increasing  the  overlap  between  bound  and  unbound  states.  For  ot  and 

+ ^ 

NO  more  definitive  data  are  available,  negating  the  need  for  this  assumption. 

Dissociative  recombination  coefficient  El  for  CO*  is  based  on  results 
at  300°K  cited  in  a review  by  Biondi  (Ref.  12)  and  a T“*  ion  temperature  de- 
pendence. 

Coefficient  E2  for  O*  is  based  on  the  experimental  work  of  Kasner  and 
Biondi  (Ref.  13). 

Coefficient  E3  for  the  process  NO*  + e — ► N + O has  been  obtained  from 

the  results  at  300  K cited  by  Biondi  (Ref.  12)  combined  with  high  temperature 

-0  5 

results  cited  by  Hansen  (Ref.  14),  influenced  by  the  T * ion  temperature  de- 
pendence for  coefficient  E2.  For  the  process  NO*  + e —►NO,  actually  corre- 
sponding to  hydrated  ion  recombination,  we  have  employed  the  value  1.0  x 10  ^ 

-0  5 - 1 

u^  ' suggested  by  the  work  of  Mitra  and  Rowe  (Ref.  4),  refit  to  a simple  T 

ion  temperature  dependence. 

Reactions  E4  through  E21;  For  the  most  part,  the  ionic  mutual  neutrali- 
zation rate  coefficients,  their  temperature  dependence,  and  the  distribution  of 
neutral  products  remain  unknown.  Recent  merged  beam  results  at  300  K are 
however,  available  (Ref.  15)  for  coefficients  E8,  E10,  El  1,  E12,  E14  and  E15 
involving  O*  and  NO*  neutralization  with  O^.NO^  and  NO^,  and  have  been 
incorporated  into  Table  1.  The  remaining  coefficients  have  of  necessity  been 
roughly  estimated  by  scaling  of  these  few  measurements  by  the  reciprocal 
square  root  of  the  reduced  mass  of  the  interacting  ions  within  the  positive  ion 
groupings  O^,  NO+  and  CO*.  On  the  basis  of  the  empirical  observation  that 
the  mean  ratio  of  the  measured  O*  mutual  neutralization  coefficients  to  the 

.j-  ^ 

corresponding  NO  coefficients  is  on  the  order  of  the  corresponding  ratio  of 
the  ion-free  electron  recombination  coefficients  at  300  K,  the  completely  un- 
known ion-ion  rates  involving  CO*  were  estimated  on  the  basis  of  the 
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corresponding  ratio  of  the  CO*  to  O*  dissociative  recombination  coefficients, 
also  scaling  via  the  reduced  mass  effect. 

As  a final  caveat,  it  is  necessary  to  point  out  that  three-body  Thompson 
ion-ion  recombination  (i.e.,  A*  + B"  + M A + B + M)  has  not  yet  been  incor- 
porated into  Table  1 but  may  play  a significant  role  in  the  ion -ion  neutraliza- 
tion kinetics  at  the  relatively  high  pressure  of  interest  in  practical  high  energy 
discharge  lasers.  Theoretical  work  by  Flannery  (Ref.  16)  indicates  that  at 
pressures  on  the  order  of  an  atmosphere  ternary  recombination  may  be  more 
rapid  than  binary  mutual  neutralization  processes. 

Excited  Species  (Table  2):  Reactions  A1  and  A4.  The  effective  rate 
coefficients  for  these  important  vibrational  pumping  reactions  are  improved 
numerical  fits  to  the  results  computed  by  Nighan  (Ref.  17)  largely  using  elec- 
tron energy  distributions  calculated  for  a 1-2-3  CC>2-N2-He  mixture.  The 
previous  numerical  fit  (i.e.,  that  given  in  Refs.  1 and  2)  for  rate  coefficient  A1 
was  not  well-behaved  at  values  of  below  ~0.45  eV,  while  that  for  A4  was 
needlessly  complex.  We  stress  again  that  the  coefficients  A1  through  A4  are 
somewhat  mixture-dependent.  For  purposes  of  exploring  the  plasma  chem- 
istry the  relatively  slight  effects  of  neglecting  this  dependence  have  been 
ignored. 

2.2.2  Combustion  Gas  Mixtures  (Hydrogenic) 

With  respect  to  C02  - CO-N2~He  laser  mixtures  with  a low  level  of 
air  contaminant  the  primary  inlet  composition  changes  in  dried  laser  mixtures 
generated  by  combustion  of  hydrocarbon  lie  simply  in  the  presence  of  a few 
percent  of  H2>  a few  hundred  parts  per  million  of  both  NO  and  unburned  hydro- 
carbon and  the  absence  of  He  (Ref.  18).  As  an  adequate  simplification  we 
assume  the  unburned  hydrocarbon  to  be  CH20  (i.e.,  formaldehyde),  a common 
combustion  intermediate. 

As  we  will  show  in  Section  3,  these  simple  changes  are  predicted  to  be 
sufficient  to  ensure  that  detrimental  levels  of  strongly  electrophilic  species 
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(primarily  HNO^  and  HNO^)  are  formed  in  the  post  electrical  pulse  dwell 
time.  The  extensions  to  the  kinetics  inputs  utilized  in  this  analysis  are  as 
follows: 
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Ionized  Species  (Table 
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• Attachment  Reactions 


Reaction  B 10;  Nighan  and  Wiegand  (Ref.  5)  have  theoretically  evaluated 
the  effective  water  dissociative  electron  attachment  coefficient,  BIO,  for  the 
products  OH  + H . To  avoid  the  necessity  of  accounting  for  H-  we  have  here 
simply  written  the  products  as  OH  + H.  Coefficient  BIO  is  a numerical  fit  to 
Nighan  and  Wiegand's  results. 


Reactions  Bll,  B12  and  B13;  Fehsenfeld,  Howard  and  Schmeltekopf 
(Ref.  19)  have  recently  reported  that  the  attachment  coefficient  of  thermal 
electrons  to  nitric  acid  at  298  K is  (5  + 3)  x 10"8  cm3  sec"1.  We  have 
accepted  this  value  and  estimated  that  the  analogous  processes  for  attach- 
ment to  HNO^  and  N^O,.  are  equally  rapid,  especially  for  relatively  hot  (i.e., 
ug  > 0.5  eV)  electrons. 

• Detachment  Reactions 

Extension  to  hydrogenic  systems  is  straightforward  via  Reactions  C9 
through  C12.  The  extensive  work  of  Ferguson,  Fehsenfeld  and  co-workers 
at  the  Boulder  Aeronomy  Laboratory  is  the  mother  lode  for  the  rate  coeffi- 
cients of  these  reactions:  C9  - Ref.  7;  CIO  and  Cl  1 - Ref.  20;  C12  - Ref.  21. 

• Charge  Transfer,  Interchange  and  Clustering  Reactions 

Similarly,  results  of  the  Aeronomy  Laboratory  are  largely  adequate  to 
extend  these  reactions  to  hydrogenic  systems,  the  sources  of  the  rate  coeffi- 
cients being  as  follows:  D36,  D37  and  D43  - Ref.  8 (coefficient  D38  has  been 
estimated  to  be  approximately  equal  to  that  for  the  generically  similar  Col) 
metathesis  D37,  leading  also  to  HCO+);  D44,  D47,  D50  and  D52  - Ref.  19; 
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Table  1 

IONIZED  SPECIES  REACTIONS?*  AND  RATE  COEFFICIENTS** 
FOR  THE  COz  ELECTRIC  DISCHARGE  LASER 


A.  Ionization  Reactions 

Rate  Coefficient 

A 1 

C02(000)  + eEB  — CO+2  + e + eEB 

6.9  x 

10'7 

„-2.0 

eeb 

A2 

c°(°)  + eEB  — CO+  + e + eEB 

4.0  x 

io-7 

tf-0.2 

eeb 

A3 

n2(°)  + eEB  -N2  + e + eEB 

3.7  x 

IQ'7 

~-0.2 

eeb 

A4 

He  + eEB  — ^ He+  + e + eEB 

6.5  x 

10'8 

p^-0.2 

eeb 

A5 

NO  + e — NO4  + e + e 

4.3  x 

10"6 

u’e4’8  exP(-15/ue) 

A6 

C02(000)  + e — *»  CO4  + e + e 

1.7  x 

10'5 

u’e6'3  exP(-19/ue) 

A7 

02  + e *•  O4  + e + e 

1 .5  x 

10'5 

ue6’3  exP(-19/ue) 

A8 

N2(0)  + e — *»N2  e + e 

7.2  x 

10'4 

u-e10  exp(-25.4/ue) 

A9 

CO(O)  + e — ► CO+  + e + e 

1.5  x 

10"5 

u"e6,°  exp(-20.3/ue) 

A10 

°2  + eEB  •*  °2  + e + eEB 

4.2  x 

10'7 

rr-  0.2 

eeb 

B.  Attachment  Reactions 

B 1 

NOz  + e — ► O"  + NO 

3.4  x 

10’10  u"e2’4  exp(-3.0/ue) 

B2 

NOz  + e — ► NO'2 

4.0  x 

10'11 

B3 

02  + e — ► O'  + O 

2.8  x 

10"6 

u‘e7-9  exp(-13.5/ue) 

B4 

03  + e — ► O"  + 02 

4.0  x 

10'11 

B 5 

NO  + e — *-  O"  + N 

3.4  x 

10'6 

u"e8,2  exp(-15.2/ue) 

B6 

C02(000)  + e — ► O'  r CO(O) 

8.4  x 

10'8 

u^8*4  exp(-11.6/ue) 

All  forward  reactions  only 

sjesjc  3 

Units  are  as  follows:  Rate  coefficients,  cm  -particle-sec  units;  Ern,  keV; 
o 1 o 1 

u_,  eV;  T,  K;  R,  cal-mole”  - K~  . (u  is  the  reduced  average  electron 
6 0 — —■ ™ ““ 

energy,  i.e.,  ug  = k Tfi). 
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1 

1 

I 

I 

B7 

O + e + M j — ► O'  + Mj 

1.0  x 

IQ'31 

B8 

°2  + e + M2  — 0'2  + M2 

1.0  x 

IQ'31 

• » 

B9 

02+e+M3  — 0-2+M3 

4.2  x 

10-27  t'1  exp(-  1200/RT) 

p 

BIO 

HzO  + e — ► OH!"  + H 

2.8  x 

10”^  a'7’2  exp(-12.9/ u ) 

I- 

B 1 1 

HNOs  + e — ► NO"2  + OH 

5.0  x 

.«-• 

B 12 

HNOz  + e — ► NO'2  + H 

5.0  x 

10'8 

a 

i 

P 

B 1 3 

N2Os  + e — ►NO'  + NOz 

5.0  x 

10'8 

li 

C.  Detachment  Reactions 

Cl 

O"  + CO(O)  — ► e + C02(020) 

5.8  x 

10'9  T'0,4 

C2 

O"  + O — ► e + 02 

1.9  x 

10-10 

C3 

O'  + NO  — ► e + NOz 

1 .8  x 

10'8T'0-75 

C4 

O'  + N — ► e + NO 

2.2  x 

10'10 

C5 

0'2  + O — ► e + 03 

3.3  x 

10-10 

C6 

0"2  + N — e + NOz 

4.0  x 

10'10 

C7  . 

03  + O — ► e + 02  + 02 

1 .0  x 

10'10 

! j 

C8 

03  + 03  ► e + 03  03 

1.0  x 

10'10 

i 

C9 

O'  + H2  — ► HzO  + e 

2.5  x 

10'9  t-0.25 

i-i 

CIO 

OH"  + H — ► HzO  + e 

1 .8  x 

10'9 

l 

i 

Cll 

OH!"  + O — ►-  HOz  + e 

2.0  x 

lO'10 

I 

u 

C12 

0'2  + H — ► HOz  + e 

1.5  x 

10'9 

1 

D.  Charge  Transfer,  Interchange 

and  Clustering  Reactions 

I 

D 1 

He+  + C02(000)  — ► CO+  + O + He 

1 .2  x 

10'9 

D2 

n2  + co2(ooo)  — ► co2  + N2(0) 

1 .0  x 

10'9 

1 n 

H 

D3 

CO+  + C02(000)  — ► C02  + CO(O) 

1 .2  x 

10"9 

LJ 

n 

D4 

C02  + Oz  — ► 0+  + C02(000) 

1.1  X 

10' 14  T exp(1800/RT) 

1 i 

■ u 

f* 

15 

■ * 

A 



L 

■1 

1.2  x 10 


I 


D6 

CO*  + O — ► O*  + CO(O) 

2.6 

X 

10 

D7 

0+  + NO  — ► NO+  + Oz 

4.5 

X 

10 

D8 

0+  + N — ► NO+  + O 

1.8 

X 

10 

D9 

o'  + o3  — ► o'3  + O 

5.3 

X 

10 

DIO 

o'  + no2  — ► no'2  + O 

1.2 

X 

10 

Dll 

O’  + Oz  + — ► 0'3  + Mj 

3.0 

X 

10 

D12 

O'  + NO  + — ► NO'2  + 

3.0 

X 

10 

D13 

O'  + C02(000)  + — ► CO"3  + Mj 

5.0 

X 

10 

D 14 

O'  + O,  + M . — ► O'  + M. 
2 4 3 4 

2.0 

X 

10 

D15 

°2  + °3  ~ °3  + °2 

3.0 

X 

10 

D 16 

°2  + N°2  N°2  + °2 

1.2 

X 

10 

D17 

0'2  + O + M L — - 0"3  + 

3.0 

X 

10 

D18 

0'2  + NO  + Mx  — ► NO'3  + M j 

3.0 

X 

10 

D19 

o'+o  — o'+o2 

1.0 

X 

10 

D20 

0'3  + C02(000)  — •-  CO"3  + o2 

6.0 

X 

10 

D21 

0'3  + NO  — ► NO'3  + O 

1.0 

X 

10 

D22 

°3  + N°2  N°3  + °2 

2.8 

X 

10 

..  D23 

CO'3  + O — ► 0"2  + CO2(020) 

8.0 

X 

10 

D24 

CO'3  + NO  — ► NO'2  + C02(020) 

1.8 

X 

10 

D25 

CO'3  + NOz  — ► NO'3  + C02(000) 

2.0 

X 

10 

D26 

NO'2  + no3  — ► NO'3  + no2 

3.0 

X 

10 

D27 

N°2  + °3  — N°3  + °2 

1.8 

X 

10 

D28 

O-  + O-  + M.  O.  + M. 

2.0 

X 

10 

I 


I 


V 


16 


10 

10 

10 

10 

10 

9 

28 


■28  T-1 


•25  T-1 


• 30 

• 10 
•9 


•28  T-1 


•28  T-1 


10 

10 

11 

10 

11 

11 

10 

10 

11 


■27  t-1.5 


D29 

o"2  + C02(000)  + m3  — ► co4  + m3 

5.2 

X 

10-29 

D30 

04  + 0 — o"3  + o2 

4.0 

X 

10"10 

D31 

°4  + c02(00°)  — *"  c°4  + °2 

4.8 

X 

10"10 

D32 

+ NO  — ► NO'  + Oz 

2.5 

X 

10"10 

D33 

co"  + 0 — ► co"  + o2 

1.5 

X 

10-10 

D34 

CO"4  + NO  — ► NO"3  + C02(000) 

4.8 

X 

10"11 

D35 

CO'4  + 03  — ► 0"3  + C02(000)  + o2 

1.3 

X 

10"l° 

D36 

C02  + H — ► HCO+  + O 

6.0 

X 

10"10 

D37 

C02  + H2  — ► HCO+  + H 

1.4 

X 

10"9 

D38 

C02  + H20  — ► HCot,  + OH 

1.4 

X 

10"9 

D39 

HCO+  + NO  — ► NO+  + CHO 

1.2 

X 

10"10 

D40 

HCO+  + HzO  — ► H30+  + CO(l) 

3.0 

X 

10"9 

D41 

HCO+  + CO(O)  — ► HCO+  + C02(020) 

3.0 

X 

10"9 

D42 

HCO+  + HzO  — »-H30+  + CO2(020) 

3.0 

X 

10"9 

D43 

O"  + HzO  — ► OH"  + OH 

1.4 

X 

10"9 

D44 

O"  + HN03  — *-  NO"3  + OH 

3.0 

X 

0 

1 

0 

►H 

D45 

OH"  + NOz  — rNO'2  + OH 

1.9 

X 

10"9 

D46 

02  + OH  — *■  OH"  + Oz 

6.0 

X 

10-10 

D47 

02  + HN03  — ► NO"  + HOz 

2.8 

X 

10"9 

D48 

O'  + H — •-  OH"  + Oz 

8.4 

X 

10-10 

D4  9 

NO"2  + H — •-  OH"  + NO 

3.7 

X 

10"10 

D50 

NO"  + HN03  — ► NO“3  + HNOz 

1.6 

X 

10"9 

D51 

CO"3  + H — ►OH"  +CO2(000) 

1.7 

X 

10"10 

D52 

CO"3  + HN03  — ► NO"3  + OH  + C02(000) 

00 

• 

0 

X 

10"10 

D53 

CO"  + H — ► CO"3  + OH 

2.2 

X 

10"10 
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E.  Recombination  Reactions 


El 

Cot;  + e — ► CO(l)  + O 

2.0 

X 

10’3  T_1‘° 

E2 

0+  +e  — ► O + O 

6.0 

X 

10"7  t“0,5 

* 

E3 

NO+  + e — ► N + O 

1.2 

X 

10"6  T-0,5 

E3 

NO+  + e — ► NO 

3.0 

X 

l0-4  t-1.0 

E4 

co2  + o'2  — ► co(i)  + o2  + O 

6.0 

X 

o 

1 

E5 

C02  +CO"3 — ► C02(020)  +O+CO2<020) 

5.0 

X 

r- 

i 

o 

rH 

E6 

co2  + no"2  — ♦ co(i)  + o + no2 

6.0 

X 

10'7 

E7 

C02  + NO"3  — ► CO(l)  + o + no3 

5.0 

X 

O 

l 

-vj 

E8 

°2  + 02-^°2+0+0 

4.0 

X 

10'7 

E9 

0+2  + CO"3  — ♦ 02  + O + C02(020) 

3.0 

X 

*—• 

O 

l 

E10 

o+  + NO'2  — •-  0+0+  no2 

4.0 

X 

10"7 

Ell 

o+2  + NO"3  — ► o + o + no3 

1.3 

X 

10"7 

E12 

NO+  + 0"2  — ► N + O + 02 

6.0 

X 

r^- 

i 

o 

E 1 3 

NO+  + CO"3  —♦NO  + O + CO2(020) 

6.0 

X 

o 

1 

E 14 

NO+  + NO'2  — ► N + O + N02 

5.0 

X 

10"7 

El  5 

NO+  + NO'3  — ► N + O + N03 

00 

b 

X 

»— * 
o 

1 

E16 

C02  + O^  — ♦ C02(020)  + 02  + Oz 

o 

• 

LD 

X 

10"7 

E17 

Cot,  + CO^  — ♦C02(020)  + C02(020) +02 

5.0 

X 

10"7 

E 1 8 

°2  t °4  °2  + °2  + °2 

3.0 

X 

10'7 

E19 

02  + CO^  — 02  + C02(020)  + 02 

3.0 

X 

10'7 

E20 

NO+  + 0‘4  — ♦ NO  + 02  + 02 

6.0 

X 

10'7 

E21 

NO+  + CO'4  —♦NO  + C02(020)  + 02 

6.0 

X 

10'7 

u 


u 


u 


li 


*s|e 


Dry 

n 

Moist 
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E22 

HCO+  + e — ► H + CO(l) 

1.0  x 

IO-5 

T- 1 u~0*^ 
e 

E23 

HCO+  + e — *■  H + C02(020) 

4.0  x 

10‘5 

T'1  u"0,5 
e 

E24 

H30+  + e — r-  H + H + OH 

4.0  x 

io-6 

T-0.5  u-0 
e 

E25 

HCO+  + 0'2  — •-  CO(l)  + H + 02 

3.0  x 

IO-7 

E26 

HCO+  + CO’  — ► CO(l)  +CO2(020)  +OH 

3.0  x 

10"7 

E27 

HCO+  + CO^  — ► CO(l)  +CO2(020)  +HOz 

3.0  x 

10-7 

E28 

hco2  + o2  — co2(020)  + h + o2 

1 .2  x 

10-6 

E29 

HC02  + C0‘3 — ► C02(020)  +CO2(020)  + OH 

1 .0  x 

10 -6 

E30 

HC02  + CO^  — ► C02(020)  + C02(020)  + HOz 

1.0  x 

10 ‘6 

E31 

h3o+  + o2  — h2o  + h + o2 

2.4  x 

10 -6 

E32 

H30+  + CO"3  — ► HzO  + C02(020)  + OH 

2.4  x 

10-6 

E33 

HsO+  + CO'4  — ► HzO  + C02(020)  + HOz 

2.4  x 

10-6 

E34 

H30+  + NO'3  — ► H20  + N02  + OH 

2.4  x 

10-6 

E35 

HsO+  + NO'2  — ► H20  + N02  + H 

2.4  x 

io-6 

M, 


M. 


M. 


Catalytic  Species 

All  species 

All  species  except  O,  O 2 and  C>3 


o(o2,o3 


= CL 
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D45  - Ref.  20;  D48,  D49,  D51  and  D53  - Ref.  21.  The  flame-ion  mass  spec- 
trometry results  of  Hurle,  Sugden  and  Nutt  (Ref.  22)  are  the  source  of  coeffi- 
cient D39.  The  rate  coefficients  for  the  proton  transfer  reactions  D41  and 
D42  have  been  estimated  to  be  approximately  equal  to  that  for  the  analogous 
reaction  D40,  the  rate  coefficient  for  which  is  that  from  the  drift  tube  study 
by  Pritchard  and  Harrison  (Ref.  23).  Finally,  coefficient  D46  for  the  charge 
transfer  reaction  between  and  OH  has  been  approximated  as  the  geometric 
mean  of  the  coefficients  for  the  analogous  charge  transfer  reactions  of  OT, 
with  03  and  NO^,  i.e.,  D15  and  D16. 

• Recombination  Reactions 

The  dissociative  ion-free  electron  recombination  coefficients  E22,  E23 
and  E24,  for  HCO  , HCO^  and  H^O*,  respectively  are  based  on  the  room  tem- 
perature results  cited  in  the  survey  by  Biondi  (Ref.  12)  of  the  effects  of  ion 
complexity  on  recombination  and  the  following:  Coefficient  E22,  for  HCO+, 
has  been  assigned  the  same  ion  temperature  dependence  as  coefficient  El, 
for  CO^;  Reaction  E23,  for  HCO*,  has  been  estimated  simply  to  be  approxi- 
mately twice  as  rapid  as  the  corresponding  dissociative  recombination  reaction 
El  involving  the  less  complex  CO*  ion;  Flame  temperature  results  by  Calcote 
et  al  (Ref.  24)  have  been  combined  with  the  room  temperature  result  cited  by 
Biondi  (i.e.,  ~1  x 10  ^)  to  infer  a ~T  ion  temperature  dependence,  which 
together  with  the  flame  results  and  the  generically  assumed  u”®'^  dependence 
leads  to  coefficient  E24. 

Ion-ion  neutralization  coefficients  E25  through  E35  have  all  been  esti- 
mated via  the  same  empirical  procedure  as  developed  (see  Section  2.2.1)  for 
the  generically  similar  coefficients  E4  through  E21,  viz.,  scaling  by  the  re- 
duced mass  effect  of  the  collision  partners  and  the  room  temperature  assess- 
ment of  the  ion-free  electron  recombination  coefficient  of  the  appropriate 
positive  ion  vis  a vis  that  for  Not 

Excited  Species  (Table  2);  Direct  kinetic  effects  on  excited  species 
of  the  addition  of  hydrogen  are  readily  and  simply  incorporated  via  inclusion 


20 


» 

I 


' 


Table  2 

EXCITED  SPECIES  REACTIONS  AND  RATE  COEFFICIENTS3 
FOR  THE  COz  ELECTRIC  DISCHARGE  LASER 


A. 

Vibrational  Pumping  Reactions 

A1 

N2(0)  + e — ► N2(l)  + e 

1.6  x 10  ^ exp(9.1  ug) 

1 + 7.0  x 10"4  exp(9. 1 u ) 

A2 

CO(O)  + e — ► CO(l)  + e 

2.4  x 

10  ^ u9'8  exp(- 1 .3l/ue’ 

1 + 5.33  ue  exp(-1.31/u  ) 

A3 

C02(000)  + e — ^ CO2(001)  + e 

7.9  x 

10~9  exp(-0.38/ue) 

A4 

C02(000)  + e — ► CO2(010)  + e 

2.0  x 

10'8  u 

e 

B. 

Electronic  Pumping  and  De-Excitation  Reactions0 

B 1 

CO(O)  + e — •-  CO*  + e 

5.3  x 

- 9 0 5 

10  7 ue*  exp(-2.6/ug) 

B2 

CO*  + Mj  — *-CO(0)+M 

1 .0  x 

lo'11" 

B3 

CO*  + C02(000)  — ► CO(O)  + O + CO(0) 

1.0  x 

O 

r— 4 
t 

o 

B4 

N2(0)  + e ►N*  + e 

1.3  x 

10"8  u~4 exp(-11.5/u 

B5 

N*  + Mj  — ^N2(0)  + Mj 

1.0  x 

io'11" 

B6 

N*  + CO2(000)  — *-N2(0)  + O + CO(O) 

1.0  x 

io-11 

B 7 

N*  + CO2(000)  — N + NO  + CO(O) 

1.0  x 

10-10 

B8 

C02(000)  + e — •-  CO*  + e 

2.5  x 

10'9  u*’6  exp(-5.1/ue) 

B 9 

CO*  + — ► C02(000)  + 

1.0  x 

io'11 

BIO 

CO*  + C02(000)  — •-  C02(000)  +0  +CO(0) 

1 .0  x 

io'10 

f f 

it 

it 

I 


aUnits  are  as  follows:  Rate  coefficients,  cm^-particle-sec  units;  u . eV; 

T,  °K;  R,  cal-mole’1  °K_1.  6 

^All  reversible  reactions;  kr  = kf/K.  For  the  vibrational  pumping  reactions, 
K is  the  equilibrium  coefficient  for  the  essentially  thermo -neutral  forward 
reaction  involving  hot  electrons. 

Q 

All  forward  reactions  only. 


y 
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C.  Vibrational  Relaxation  Reactions 


Cl  CO2(110)  + C02(000)  — ►CO2(100)  + CO2(010) 
C2  CO2(030)  + C02(000)  — ► CO2(100)  + CO2(010) 
C3  C02(030)  + CC>2(000)  — »-C02(020)  + CO2(010) 
C4  CO2(100)  + C02(000)  — ► CO2(0 10)  + CO2(010) 
C5  C02(020)  + C02(000)  — ►CO2(010)  + CO2(010) 
C6  C02(001)+M5 — ^CO2(110)+M5 
Cl  CO2(001)  + Mg  — ► CO2(110)  + Mg 
C8  CO2(001)+M5 — »- C02(030)  + M5 
C9  CO2(001)  + Mg  — «-  C02(030)  + Mg 
CIO  CO2(110)  + Mg  — ► C02(030)  + Mg 
Cll  CO2(110)  + Mg — ►C02(020)+Mg 
C 12  CO2(110)  + Mg  — ► C02(020)  + Mg 

C 1 3 CO2(110)  + Mg — ►CO2(100)+Mg 

C 14  C02(030)  + Mg — ►C02(020)+Mg 

C 15  C02(030)  + Mg — ►CO2(100)+Mg 

C 16  CO2(100)  + Mg — ►C02(020)+Mg 

C 1 7 CO2(100)  + M?  — *-  CO2(010)  + M? 

C 18  CO2(100)  + He — ►CO2(010)+He 

C19  CO2(100)  + CO(O)  — CO2(010)  + CO(O) 

C20  CO2(020)  + M?  — ► CO2(010)  + M? 

C21  CO2(020)  + He — •- CC>2(0  10)  + He 

C22  C02(020)  + CO(O)  — CO2(101)  + CO(O) 


1.25  x 10‘13  T0,5 
1.8  x 10"15  T0,5 
3.1  x 10'13  T0,5 


4.0  x 10 


-13 


1.4  x 10-12  T°*5 

1.1  x 10"27  T4,8  exp(1484/RT) 
1.9  x 10"31  T5,8  exp(2436/RT) 

8.1  x 10‘31  T5,6  exp(1484/RT) 

1.4  x 10"34  T6,6  exp(2436/RT) 

4.3  x 10  Ta#d 

4.5  x 10'27  T4*2  exp(-903/RT) 

8.8  x 10"20  T2-5  exp(-4410  /RT) 

8.6  x 10"24  T3*8  exp(-549/RT) 

9.3  x 10"22  T3,3  exp(- 1230  RT) 

1.1  x 10'21  T3,0  exp(-  1060  RT) 

8.0  x 10‘18  T1*5 

5.6  x 10"22  T3,3  exp(-  1480/RT) 

1.8  x 10'21  T3,0  exp(843/RT) 

8.4  x 10'10  T"1,0 

2.1  x 10"21  T3,2  exp(-  1350/RT) 

3.8  x 10'21  T3,0  exp(843/RT) 

1 -7  10-9  rp“  1 .0 

1.7x10  T 


a 


All  reversible  reactions. 
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C23 

C02(010)  + M?  — •-  C02(000)  + M? 

3.4  x 

IQ'26 

T4,2  exp(l  1 30/RT) 

C24 

CO2(010)  + He  — *-  C02(000)  + He 

9.9  x 

io-22 

T3,0  exp(843/RT) 

C25 

CO2(010)  + CO(O)  — *-C02(000)  + CO(O) 

4.5  x 

io-10 

o 

• 

1 

H 

C26 

C02(001)  + N2(0)  — C02(000)  + N2(l) 

8.3  x 

io'12 

t-0.5 

C27 

C02(001)  + CO(O)  — ► CO^tOOO)  + CO(l) 

1.4  x 

io"17 

rpl  .65 

C28 

N2(l)  + CO(O)  — *-N2(0)  + CO(l) 

2.7  x 

io-18 

t1.5 

C29 

CO2(010)  + H20  — ► C02(000)  + h2o 

3.2  x 

10-13 

exp(22.9/T1/3) 

C30 

CO2(001)  + HzO  — «-  C02(030)  + h2o 

4.0  x 

IO-13 

C31 

N2(l)  + HzO  — ►N2(0)  + h2o 

1.1  X 

io-10 

exp(-68.9/T1/3) 

C32 

CO(l)  + HzO  — ► CO(O)  + h2o 

3.1  x 

io-10 

exp(-65.0/T1/3) 

C33 

CO2(010)  + H2  — C02(000)  + H2 

3.0  x 

io-12 

C34 

CO2(001)  + H2  — ^CO2(110)  + H2 

1.7  x 

10"17 

T 1 .5 

Catalytic  Species 

= All  species 

M,-  = All  species  except  He  and  N2 

= 1.5  He;  All  others;  1.0 

= 0.5  N2;  CO  = He  = 0;  All  others;  1.0 

] 


of  the  dominant  relaxation  rates  of  the  upper  and  lower  laser  levels  of  CO^ 
by  Reactions  C33  and  C34.  Rate  coefficients  C33  and  C34  are  the  result 
of  curve  fits  to  the  recent  data  survey  and  theoretical  evaluations  by  Strieker 
(Ref.  25). 

Free  Radical  Species  (Table  3):  For  the  most  part  the  additional  re- 
actions and  rate  coefficients  required  to  extend  the  CC^-CO-^-He  free 
radical  species  hinetics  given  in  Ref.  2 to  combustion  gas  mixtures  (i.e., 
Reactions  A8  through  A17  and  B 14  through  B36)  have  been  taken  directly 
from  the  current  Lockheed  LAMP  reaction  models  used  fairly  extensively 

j|c 

for  analysis  of  chemical  lasers  and  rocket  exhaust  plumes.  The  few  exceptions 
are  discussed  below. 


• Ternary  Recombination  Reactions.  The  ternary  rate  coefficient  for 
Reaction  A 14  leading  to  HNO3  is  that  adopted  in  the  theoretical 
stratospheric  model  developed  by  Shimazaki  and  Ogawa  (Ref.  27). 
Reactions  A15,  A16  and  A17  are  normally  ternary  in  stratospheric 
applications  (see,  e.g.,  Ref. 28).  At  pressures  on  the  order  of  atmos- 
pheric, however,  they  follow  essentially  binary  behavior,  being  at 
their  high  pressure  limit  (Ref. 29).  Coefficients  A15.A16  and  A17 
are  the  recommended  survey  values  given  by  Baulch  et  al.,  (Ref. 29) 
for  these  rate  coefficients  in  their  binary  limit. 

• Binary  Reactions.  Rate  coefficient  B19  for  the  H+NO2  reaction  is 
the  recommended  survey  value  given  by  Baulch  et  al.,  (Ref.29). 
Coefficients  B 17,  B22,  B23,  B24,  B27  and  B30  are  those  adopted  by 
Wuebbles  and  Chang  (Ref.  30)  in  their  study  of  the  temporal  behavior 
of  the  stratospheric  distribution  of  miner  chemical  constituents. 
Coefficient  B16  for  the  reaction  between  OH  and  HNO2  has  been 
estimated  as  approximately  equal  to  that  for  the  analogous  reaction, 
B17,  between  OH  and  HNO^. 


These  models  are  kept  as  current  as  possible  by  incorporation  of  recent 
results  as  these  are  made  available,  the  most  recent  documentation  being 
Ref. 26.  Validation  via  correlation  with  data  has  largely  been  successful. 
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Table  3 


FREE  RADICAL  SPECIES  REACTIONS'1  AND  RATE  COEFFICIENTS6 
FOR  THE  C02  ELECTRIC  DISCHARGE  LASER 


A. 

Ternary  Recombination  Reactions 

A1 

N + 

N + M9— N2(l)  +M9 

A2 

N + 

O + Mg  ►NO  + Mg 

A3 

O + 

CO(0)  + M1q  — C02(020)  + 

A4 

O + 

O + Mg—  02  + Mg 

A5 

O + 

°2  + M9  — *■  °3  + M9 

A6 

o + 

NO  +Mg — ►N02  +Mg 

A7 

o + 

NOz  + Mg  ►N03  + Mg 

A8 

H + 

H + Mg  — ► H2  + Mg 

A9 

O + 

H + Mg  — ► OH  + Mg 

A10 

H + 

OH  + Mg  — ► HzO  + Mg 

All 

H + 

02  +Mg  "“H02  + Mg 

A12 

CHO  + — ► CO(0)  + H + 

A13 

ch2 

,0  + Mg  — ► CO(0)  + H2  + Mg 

A14 

h°2 

, + NO  + Mg  — ► HN03  + M 

A15 

OH 

+ NOz  — ► hno3 

A16 

OH 

+ NO  — ► HN02 

A17 

NO. 

, +N°3— N2°5 

8.3  x 

1.8  x 

2.0  x 

3.8  x 

1.0  x 

4.0  x 

1 .8  x 

2.8  x 

1.0  x 

6 . 1 x 

2.1  x 

1.2  x 
3.5  x 
3.5  x 

1.0  x 

2.0  x 

4.0  x 


Rate  Coefficients 


10 

10 

10 

10 

10 

10 

10 

10 

10 

10 

10 

10 

10 

10 

10 

10 

10 


exp(1000/RT) 

"31  T’0,5 

’33  exp-(4000/RT) 

*30  T-1  exp-(340/RT) 

‘34  exp(1010/RT) 

"33  exp(1880/RT) 

"32  exp(  1000 /RT) 

-30  T-1.0 

-32 

-26  t-2.0 

"32  exp(580/RT) 

"l0  exp(-  15000/RT) 

"8  exp(-35000/RT) 

-31 

*U  exp(-340/RT) 

-12 

-12 


All  reversible  reactions. 

b 3 

Units  are  as  follows:  Rate  coefficients,  cm  -particle-sec  units;  T,  °K; 

R,  cal-mole“  *-°K“ 


B.  Binary  Reactions 


B 1 

N + NO  — ►N2(0)  + O 

2.7  x 

io-11 

B2 

N + NOz  — ►NO  + NO 

1.4  x 

IQ'12 

B 3 

N + N03  — ►NO  + N02 

5.7  x 

io'13 

B4 

N + Oz  — ► NO  + O 

1.1  X 

io-14 

B 5 

N + 03  — ►NO  + Oz 

5.7  x 

io'13 

B6 

° + °3  — * °2  + °2 

1.9  x 

10-11 

B7 

O + NOz  — ► 02  + NO 

1 .6  x 

io-10 

B8 

o + no3  — ► o2  + no2 

5.7  x 

10'13 

B9 

NO  + N03  — ► NOz  + NOz 

00 

X 

io-12 

. -13 

BIO 

NO  + 03  — ► NOz  + Oz 

9.0  x 

10 

Bll 

NOz  + o3  — ► no3  + o2 

1 .2  x 

io-13 

B 12 

N02  + N03  — ► N02  + Oz  + NO 

2.3  x 

rO 

i 

O 

B 1 3 

no3  + no3  — ► no2  + o2  + no2 

5.0  x 

io'12 

B 14 

OH  + CO(O)  — ► C02(000)  + H 

1.1  x 

10"19 

B 15 

OH  + H2  — ►HzO  + H 

1 .0  x 

io-17 

B 16 

OH  + HNOz  — ► HzO  + NOz 

3.5  x 

10“13 

B 1 7 

OH  + HN03  — ► HzO  + N03 

3.5  x 

IO-13 

B 18 

OH  + N — ►NO  + H 

5.3  x 

10'11 

B 19 

H + NOz  — ► OH  + NO 

6.0  x 

io-10 

B20 

OH  + O — ► H + 02 

4.0  x 

10'11 

B21 

OH  + OH  — ► HzO  + O 

1.0  x 

10-11 

B22 

H + 03  — ► OH  + 02 

2.6  x 

io-11 

B23 

OH  + 03  — ► HOz  + Oz 

1.6  x 

10"12 

T exp-(6300/RT) 
exp-(4600/RT) 

ip-0.5 


exp-(2400/RT) 
exp-  (4900/RT) 
exp-(3200'RT) 
exp-(6000/RT) 

T2  exp(l600  RT) 
T2  exp(-2900  RT) 
exp(-600/RT) 
exp(-600/RT) 

exp(-  1500/RT) 

exp(-  1 100/RT) 

exp(-2000/RT) 
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B25 
B 26 
B27 
B28 
B29 
B30 
B31 
B 32 
B33 
B34 
B35 
B36 


-^no2+no2+o2 
H02  + H — ► OH  + OH 

ho2+h-~h2  + o2 

H02  + NO  — ► NOz  + OH 
HOz  + O — ►OH  + 02 
HOz  + OH  — ► HzO  + 02 
HOz  + 03  — ► OH  + 02  + Oz 
CHzO  + H — ► CHO  + H2 
CH20  + O — ► CHO  + OH 
CH20  + OH  — ► CHO  + H20 
CHO  + O — ► CO(0)  + OH 
CHO  + OH  — ► CO(O)  + H20 
CHO  + Oz  — ► CO(O)  + H02 


1 .0  x 
4.2  x 

4.2  x 

2.0  x 

8.0  x 

8.3  x 

1.0  x 

2.2  x 
1 .6  x 

9.0  x 

2. 1 x 
2.1  x 

8.3  x 


10 


10 


-14 


-10 


10 


10 


-11 


-13 


10 


11 


10 


10 


-11 


-13 


10 


11 


10 


•13 


10 


-13 


10 


-10 


10 


•10 


10 


11 


exp(-  1900/RT) 
exp(-700/RT) 

exp(-  1000/RT) 
exp(-  1000/RT) 
exp(-2500/RT) 
exp(-3800/RT) 


exp(-  1600/RT) 


Catalytic  Species 

Mg  = 2 C02>  He;  All  others:  1.0 
M 1Q  =3  C02,  1.5  CO;  20  Oz;  All  others;  1.0 
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Section  3 
RESULTS 


3.1  COMPARISON  WITH  EXPERIMENT 

In  order  to  put  results  and  possible  recommendations  derived  from  the 
theoretical  model  discussed  in  Section  2 into  perspective  it  was  felt  desirable 
to  first  test  and  validate  the  model  against  some  available  experimental  data. 
Suitable  results  for  comparison  were  found  in  a paper  by  Dezenberg,  Cason 
and  Huff  (Ref.  31)  which  presents  single  pulse  performance  data  for  an  electron 
beam-sustained  electric  discharge  laser.  These  data  have  the  additional  ad- 
vantage that  they  were  obtained  using  the  standard  CO^/N^/He  = 1/2/3 
mixture  which  was  previously  analyzed  by  the  present  authors  in  connection 
with  a study  of  a continuous  wave  electron  beam-sustained  laser  (Refs.  1,  2). 
The  first  theoretical  results  obtained  in  the  present  study  pertain  to  the  case 
of  5 atm  initial  cavity  pressure  as  discussed  in  Ref.  31.  The  E/p  value  of  3.5 
kV  cm'*  atm"*  as  reported  there  then  implies  an  electric  field  strength 
of  17.5  kV/cm,  and  — using  the  given  electrode  separation  of  4 cm  — an  applied 
drift  field  potential  of  70  kV.  Unsaturated  gain  measurements  were  made 
using  a nominal  post-foil  electron  beam  current  density  equal  to  220  mA/cm^, 
which  was  also  used  in  the  present  calculations.  Oscillator  performance 
calculations  with  the  present  model  considered  optical  cavity  data  as  indicated 
in  Refs.  31  and  32  using  an  8%  reflectivity  for  the  Na  Cl  output  flat  (Ref.  33). 
Theoretical  results  for  this  case  are  summarized  and  compared  with  experi- 
ment in  Figs.  1 through  4 and  Table  4.  Figure  1 shows  the  temporal  behavior 
of  the  post-foil  electron  beam  current  density  assumed  in  our  calculations, 
while  Fig . 2 shows  the  drift  field  voltage  as  observed  in  the  experiment  and  as 
used  in  our  calculations.  The  pulse  width  used  in  the  calculations  was  set  at 
1.6  Jisec  as  indicated  by  the  sharp  drop  in  the  experimentally  observed  distri- 
butions of  gun  current  (see  Ref.  31)  and  drift  field  voltage  at  that  time.  While 
the  present  calculations  assume  an  instantaneous  cutoff  in  both  electron  beam 
current  density  and  drift  field  voltage,  the  experiment  shows  a finite  rate  of 
decay  at  least  for  the  latter,  indicating  the  possibility  of  some  pumping 


Fig.  2 - Temporal  Behavior  of  Drift  Field  Voltage  at  5 atm 
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Table  4 

PERFORMANCE  COMPARISON 
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Theory 
(This  Work) 

Experiment 
(Ref.  31) 

Initial  Pressure  (atm) 

5.0 

5.0 

Cavity  Volume  (lit) 

4.0 

4.0 

Drift  Field  (kV/cm-atm) 

3.5 

3.67 

Output  Energy  Density 
( J/lit-atm) 

10.2 

12.5 

Pulse  Energy  (J/pulse) 

204 

250 

beyond  the  nominal  duration  of  the  pulse.  Further  evidence  for  this  is  seen 
in  the  temporal  behavior  of  the  resulting  drift  field  current  shown  in  Fig.  3. 
While  the  agreement  between  theory  and  experiment  is  excellent  over  the 
major  part  of  the  pulse  width,  the  calculated  drift  field  current  vanishes  at 
the  end  of  the  pulse  simply  as  a result  of  the  vanishing  drift  field  potential. 
Unsaturated  gain  distributions  as  measured  and  as  calculated  are  compared 
in  Fig.  4.  It  is  interesting  to  note  that  Dezenberg,  Cason  and  Huff  state  in 
their  paper  that  gain  coefficients  in  excess  of  5%  cm"'  were  obtained  at  5 atm, 
just  as  predicted  by  the  present  calculations.  Overall,  the  agreement  of 
measured  and  calculated  temporal  behavior  of  the  unsaturated  gain  is  con- 
sidered very  satisfactory. 

Finally,  calculated  laser  output  energy  is  compared  with  experimental 
data  as  shown  in  Table  4.  Although  the  experimental  drift  field  strength  of 
3.67  kV/cm-atm  at  otherwise  unchanged  conditions  indicates  a discharge 
voltage  slightly  higher  than  70  kV,  the  theoretical  results  are,  for  simplicity, 
based  on  the  previously  mentioned  conditions  (i.e.,  3.5  kV/cm-atm).  Both 
theoretical  output  energy  density  and  pulse  energy  are  approximately  18% 
lower  than  their  measured  counterparts.  Considering  the  functional  relation- 
ship between  power  input  and  E/p  only  (see  Eq.  (16)),  a roughly  15%  lower 
value  for  the  calculated  laser  output  would  be  expected  if  it  can  be  assumed 
that  efficiency  is  unaffected.  Considering  in  addition  the  differences  in 
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post-pulse  temporal  behavior  of  electron-beam  current  and  drift  field  voltage, 
the  difference  between  theory  and  experiment  appears  entirely  accounted  for. 

It  can  therefore  be  concluded  that  theory  and  experiment  are  basically  in 
agreement. 

3.2  SINGLE  PULSE  BEHAVIOR  OF  SEVERAL  GAS  MIXTURES 

Using  identical  electric  discharge  conditions  an  investigation  was  made 
of  the  theoretical  single  pulse  performance  of  two  mixtures  containing  hydrogen 
in  order  to  compare  with  the  standard  CO^/N^/He  = 1/2/3  mixture.  The  first 
mixture  was  of  clean  CO^/N^/H^  = 1/3/0. 1 composition  while  the  second  mixture 
was  representative  of  a laser  gas  that  could  conceivably  be  obtained  by  hydro- 
carbon combustion,  namely  CO^/CO/N^/H^  = 1/0.25/5.5/0.125  including  0.15% 
of  and  CH^O  each,  as  well  as  0.015%  NO  (Ref.  18). 

For  all  the  gases,  initial  conditions  in  the  laser  cavity  were  assumed  to 

be  200  K and  1 atm.  Post-foil  electron  beam  current  density  and  effective 

2 

beam  electron  energy  were  set  at  100  mA/cm  and  100  keV,  respectively.  The 
applied  discharge  potential  was  assumed  to  be  25  kV  across  the  5 cm  gap  be- 
tween the  discharge  electrodes,  resulting  in  an  electric  field  of  5 kV/cm  and 
an  E/N  value  of  1.36  • 10  V cm  . Pulse  widths  for  the  electron  beam  and 
the  discharge  were  8 /i sec  and  10  flsec,  respectively,  both  being  initiated  simul- 
taneously. Only  square  pulses  were  investigated.  Discharge  electrode  dimen- 
sions were  assumed  to  be  5 x 25  cm,  resulting  in  a discharge  volume  of  0.625 
liters.  The  conditions  assumed  are  representative  for  the  small  scale  circu- 
lator planned  by  the  Army  (Ref.  34). 

Calculated  drift  field  current  density  distributions  are  shown  in  Fig.  5. 
Corresponding  values  for  the  total  drift  field  energy  input  are  590,  510  and 
755  J lit-1  atm'1  for  the  three  mixtures  C02/N2/He,  C02/N2/H2  and  CO 2/ 
CO/N2/H2,  respectively.  Temporal  behavior  of  the  secondary  electron 
energy,  of  the  gas  temperature  and  the  unsaturated  gain  are  shown  in  Figs.  6, 

7 and  8,  respectively.  It  is  interesting  to  note  that  the  gain  begins  to  decay 
in  all  three  gas  mixtures  as  soon  as  the  temperature  exceeds  about  350  K. 


34 


Plasma  chemistry  effects  in  the  three  gas  mixtures  are  shown  in  Figs. 

9,  10  and  1 1 in  terms  of  temporal  distributions  of  the  major  species  being 
formed  in  the  discharge.  Figure  9 displays  results  for  the  CO^/N^/He  mix- 
ture, with  neutral  species  being  shown  in  Fig.  9a,  electrons  and  negative  ions 
in  Fig.  9b,  and  positive  ions  in  Fig.  9c.  For  the  CO^/N^/He  mixture,  the  major 
electrophilic  species  to  be  expected  are  and  NO^.  In  a single  pulse  dis- 
charge only  the  latter  is  being  formed  in  any  noticeable  amount  as  seen  in 
Fig.  9a,  and  at  this  level  NO^  apparently  has  no  appreciable  effect  on  the  free 
electron  concentration  shown  in  Fig.  9b.  The  dominant  neutral  species  formed 
in  the  discharge  are  atomic  oxygen,  carbon  monoxide,  plus  atomic  nitrogen  and 
nitric  oxide  in  equal  amounts.  A very  small  amount  of  oxygen  and  nitrogen 

dioxide  is  also  predicted.  The  dominant  ions  formed  are  CO,  and  CO*  as 

+ i c 

seen  from  Figs.  9b  and  9c.  Note  that  CO^,  its  distribution  practically  paralleling 
that  of  the  free  electrons,  almost  entirely  accounts  for  overall  charge  neutrality 
in  the  plasma. 

Plasma  chemistry  results  for  the  CO^/^/H^  mixture  are  shown  in 
Fig.  10.  As  seen  in  Fig.  10a,  atomic  hydrogen  now  is  by  far  the  dominant 
neutral  species  being  generated  in  the  discharge.  CHO,  N,  NO  and  O all 
reach  levels  of  a few  parts  per  million  at  the  end  of  the  discharge,  and  every- 
thing else  is  formed  only  in  even  smaller  amounts.  The  formation  of  NO^ 
barely  reaches  the  one  part  per  billion  level  at  the  end  of  the  discharge.  Thus, 
just  as  for  the  CO^/N^/He  mixture  discussed  previously,  the  essential  absence 
of  major  electrophilic  species  during  the  first  pulse  results  in  a steady  in- 
crease of  the  free  electron  density  and  the  drift  field  current  density  as  shown 
in  Figs.  10b  and  5,  respectively,  until  the  electron  beam  pulse  terminates  at 
8 (isec.  CO,  is  again  the  most  prominant  negative  ion,  however,  HCO*  now 
is  the  dominant  positive  ion  with  CO^  at  a level  approximately  two  orders 
of  magnitude  lower  (see  Figs.  10b  and  10c). 

It  is  important  to  realize  that,  of  course,  the  formation  of  no  species 
can  be  predicted  in  the  type  of  calculations  being  discussed  here  unless  it  is 
provided  for  in  the  set  of  reactions.  Theoretical  results  like  these  therefore 


Time  (/isec) 

Fig.  10b  - CO2/N2/H2  = 1/3/0. 1;  Predicted  Single  Pulse  Temporal  Behavior 
of  Major  Discharge-Generated  Negatively  Charged  Species 
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must  be  reassessed  constantly  as  to  their  validity.  It  is  believed  that  the 
results  shown  in  Fig.  10,  and  in  Fig.  10a  in  particular,  may  not  be  entirely 
correct.  The  reason  is  that  even  small  concentrations  of  certain  species 
can  lead  to  the  formation  of  other  species  which  may  become  very  important. 
Analysis  of  the  data  shown  in  Fig.  10a  led  us  to  believe  that  the  amount  of  OH 
formed  is  sufficient  to  lead,  in  conjunction  with  the  NO  available,  to  the  forma- 
tion of  HNO^  (nitrous  acid)  via  the  reaction  OH  + NO  — HNO^.  Similarly, 
an  even  smaller  amount  of  HO^  (hydroperoxyl)  when  combined  with  NO  will 
result  in  the  formation  of  HNO^  (nitric  acid)  via  the  reaction  HO^  + NO  — ► 
HNO^j.  The  importance  of  this  realization  lies  in  the  fact  that  both  HNO^  and 
HN03  are  strongly  electrophilic  and  may  therefore  exert  an  appreciable  in- 
fluence on  the  free  electron  distribution.  In  fact,  this  may  be  the  very  reason 
for  the  early  decay  of  the  drift  field  current  density  in  the  CO^/CO/N^/H^ 
mixture  as  shown  in  Fig.  5 since  the  analysis  of  this  particular  mixture  did 
include  the  formation  of  HN02  and  HN03.  This  is  shown  in  more  detail  in 
Figs.  11a  through  11c.  As  can  be  seen  in  Fig.  11a,  both  OH  and  HO^  reach 
the  level  of  one  part  per  million  very  early  in  the  discharge.  With  NO  already 
available  as  a result  of  hydrocarbon  combustion  this  leads  to  relatively  strong 
formation  of  HNO^  and  HN03<  The  other  major  electrophilic  species,  03  and 
NO^  also  reach  levels  of  one  part  per  million  and  almost  20  parts  per  million, 
respectively,  at  the  end  of  the  discharge. 

Ions  formed  in  the  CO^/CO/N^/H^  combustion  mixture  are  shown  in 
Figs,  lib  and  11c.  Among  the  negative  ions,  NO’^  is  the  only  one  displaying 
a positive  gradient  at  the  end  of  the  discharge  indicating  a certain  degree  of 
stability.  A previous  study  of  a closed  cycle  continuous  wave  electric  dis- 
charge laser  (Ref.  2)  indicated  that  NO,  was  the  most  stable  negative  ion. 

+ ^ + 

Fig.  11c  shows  HCO  to  be  the  dominant  positive  ion,  and  NO  to  be  the  second 

most  important  one.  Again,  previous  studies  indicate  that  NO*  is  likely  to  be 

the  most  stable  positive  ion  as  indicated  by  its  relatively  slow  decay  at  the 

end  of  the  pulse. 


With  calculations  presented  here  being  essentially  of  exploratory  nature 
it  is  difficult  to  reach  any  firm  judgment  about  the  relative  merits  of  the  laser 
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gas  mixtures  discussed.  Judging  by  the  peak  gain  which  was  predicted  it 
seems  that  none  of  the  two  mixtures  containing  can  match  the  perform- 
ance of  the  CC>2/N2/He  mixture.  However,  the  peak  gain  predicted  for  the 
C02/C0/N2/H2  mixture  is  sufficiently  high  to  justify  further  exploring  the 
concept  of  laser  gas  generation  via  hydrocarbon-air  combustion.  Much  de- 
pends on  what  happens  if  these  mixtures  are  pulsed  repetitively.  For  ex- 
ample, previous  studies  (Ref.  2)  have  indicated  that  if  a CO_/CO/N-/He 
mixture  is  recirculated  through  a continuous  discharge,  steady  state  per- 
formance can  easily  drop  to  50%  of  single  pulse  performance  even  with  the 
use  of  fresh  makeup  gas. 


Some  exploratory  calculations  of  repetitively  pulsing  the  CO_/CO/N_/He 

£ £ 

mixture  discussed  above  are  presented  next. 


3.3  REPETITIVE  PULSING  AND  RECIRCULATION 


One  of  the  goals  of  the  present  investigation  was  to  shed  some  light  on 
the  response  of  EDL  plasmas  when  subjected  to  recirculation  and  repetitive 
pulsing.  Being  of  primary  interest,  the  CC>2/CO/N2/H2  combustion-generated 
laser  gas  mixture  was  analyzed  according  to  the  description  given  in  Section  2. 
The  simulated  recirculation  loop  consisted  of  a discharge  pulse  of  10  jLlsec 
duration,  followed  by  an  adiabatic  expansion  process  of  100  JLlsec  duration, 
followed  by  an  isobaric  cooling  process  of  5msec  duration.  This  cycle  was 
repeated  five  times,  each  cycle  starting  with  the  initial  conditions  of  1 atm 
and  200  K,  and  using  the  gas  composition  essentially  as  obtained  at  the  end 
of  the  previous  cycle.  Three  variations  were  analyzed,  and  their  key  results 
are  tabulated  for  easy  comparison  in  Table  5.  For  Case  I,  the  laser  gas  was 
repetitively  pulsed  without  any  modification  between  cycles.  As  seen  in  Table  5, 
a drastic  reduction  in  peak  gain  and  discharge  energy  density  occurred  after  the 
first  cycle  with  an  attendant  shift  of  the  peak  gain  to  the  end  of  the  electron 
beam  pulse  (at  8 jLlsec).  A moderate  rise  in  both  peak  gain  and  discharge  energy 
density  commenced  with  the  fourth  cycle  indicating  that  after  five  cycles  the 


plasma  was  still  in  a transition  phase  rather  than 


in  a truly  repetitive  situation. 
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Table  5 

RECIRCULATION  AND  REPETITIVE  PULSING 
COz/CO/N2/H2  = 1/0.25/5.5/0.125  Mixture 


Case  Cycle 


max 

(%/cm) 


“max  InPut 

(pLsec)  ( J lit-* -atm"  *) 


Remarks 


758.0 

756.0 


No  Scrubbing 


Scrubbing  H20  Only 


Scrubbing  H20,  HN02.  HN03 


For  Case  II,  the  water  formed  in  the  discharge  (see  Fig.  11a)  was  re- 
moved at  the  end  of  each  cycle,  the  thought  being  that  any  heat  exchanger 
which  would  cool  the  gas  to  200  K would  certainly  freeze  out  any  moisture 
contained  in  the  gas.  As  the  results  show,  the  situation  is  qualitatively  the 
same  as  for  Case  I,  and  quantitatively  even  slightly  worse  as  far  as  reduction 
of  peak  gain  is  concerned.  For  Case  III,  not  only  water  but  also  HNC>2  and 
HNO^  formed  during  the  cycle  were  selectively  removed  at  the  end  of  each 
cycle,  with  the  result  that  now  the  peak  gain  and  its  location  in  time  remains 
practically  the  same  for  each  cycle  and  equal  to  the  peak  gain  of  the  initial 
pulse. 
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The  foregoing  results  clearly  establish  the  fact  that  the  formation  of 
both  nitrous  and  nitric  acid  play  a major  role  in  degrading  the  plasma  in  a 
closed  loop  EDL  using  hydrocarbon -air  combustion  gas  as  the  lasing  medium. 
Figures  12  and  13  serve  to  further  illustrate  the  results  obtained  and  dis- 
cussed above.  Both  figures  show  on  a distorted  time  scale  the  temporal  con- 
centrations of  key  contaminants  formed  during  the  five  cycles.  In  order  to 
better  visualize  what  happens  during  each  cycle,  the  first  half  of  each  interval 
on  the  abscissa  represents  the  discharge  of  10  /Ltsec  duration,  while  the  second 
half  of  each  interval  represents  the  remainder  of  the  cycle  of  approximately  5 
msec  duration.  Figure  12,  first  of  all,  shows  that  once  the  adiabatic  expansion 
is  over  all  species  gradients  essentially  vanish.  This  indicates  that  the  plasma 
is  virtually  chemically  frozen  through  the  remainder  of  the  cycle  rendering 
finite  rate  calculations  in  this  portion  of  the  cycle  largely  superfluous.  Second 
considering  the  temporal  behavior  of  the  HNO^  concentration,  which  is  gen- 
erally about  one  order  of  magnitude  higher  than  the  HNO^  concentration,  it 
can  be  assumed  that  it  is  mainly  HNO^  which  causes  the  deterioration  of  the 
gain.  Furthermore,  low  average  HNO^  concentrations  over  the  duration  of 
the  discharge  pulse  appear  to  correlate  well  with  high  peak  gain  values,  and 
vice  versa,  which  also  indicates  a close  relationship  between  gain  (or  electron 
concentration)  and  HNO^  concentration.  Otherwise,  Fig.  12  shows  again  that 
a truly  repetitive  or  steady  state  has  not  been  reached  within  five  cycles. 

This  is  contrasted  by  the  truly  repetitive  temporal  behavior  of  some  key 
species  concentrations  shown  in  Fig.  13  which  represents  the  case  where 
H^O,  HNO^  and  HNO^  are  removed  from  the  gas  after  each  cycle. 


r 

r. 


h 


! 


52 


Fig.  12  - CO2/CO/N2/H2  = 1/0.25/5.5/0.125;  Predicted  Temporal  Behavior 
of  Selected  Discharge -Generated  Species  Under  Closed  Cycle  and 
Retetitive  Pulsing  Conditions 
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Cycle  Number 

Fig.  13  - CO2/CO/N2/H2  = 1/0.25/5.5/0.125;  Predicted  Temporal  Behavior 
of  Selected  Discharge -Generated  Species  Under  Closed  Cycle  and 
Repetitive  Pulsing  Conditions  with  Selective  Removal  of  I^O,  HNO 
and  HNO,  after  Each  Cycle 
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Section  4 

CONCLUSIONS  AND  RECOMMENDATIONS 

A theory  has  been  formulated,  and  a computer  code  has  been  developed 
to  investigate  the  influence  of  discharge  and  combustion -produced  molecular 
species  on  the  performance  of  pulsed  electric  discharge  lasers.  As  com- 
pared to  some  previous  analyses  of  pulsed  CO^  electric  discharge  lasers  (for 
example,  Refs.  35  through  38),  the  significance  of  the  analysis  presented  here 
lies  in  the  coupling  of  the  detailed  plasma  chemistry  to  the  other  processes 
of  importance  in  electric  discharge  laser  cavities.  Performance  models  for 
pulsed  electric  discharge  lasers  by  Vlases  and  Moeny  (Ref.  35),  Kast  and  Cason 
(Ref.  36),  and  Douglas -Hamilton  and  Lowder  (Ref.  38)  consider  electron  and  ion 
concentrations  in  terms  of  global  coefficients  for  generation,  recombination 
and  attachment  but  otherwise  ignore  plasma  chemistry  effects.  According  to 
Wiegand  and  Nighan  (Ref.  39),  the  model  of  plasma  chemical  phenomena  de- 
veloped by  Wiegand  (Ref.  37  which,  however,  does  not  present  any  details) 
considers  approximately  300  reactions  and  follows  the  development  of  nearly 
40  neutral  and  charged  species.  However,  apparently  neither  the  effects  of 
gas  temperature  variations  nor  the  process  of  lasing  itself  was  included. 

Theoretical  predictions  obtained  from  the  present  analysis  were  com- 
pared with  experimental  data  for  a pulsed  electric  discharge  laser  using 
the  standard  CC^/^/He  = 1/2/3  mixture.  Encouraging  agreement  between 
theory  and  experiment  was  obtained  for  the  peak  gain  and  its  location  in  time, 
as  well  as  for  the  pulse  energy. 

Having  verified  the  validity  of  the  theory,  the  code  was  applied  to  a^tom- 
parative  evaluation  of  two  other  potential  laser  gas  mixtures.  Both  mixtures 
contained  H2  in  place  of  He,  and  the  composition  of  one  mixture  was  modeled 
to  closely  represent  a hydrocarbon-air  combustion -produced  laser  gas.  Con- 
sidering the  peak  unsaturated  gain  coefficient  as  a measure  of  comparison, 
both  mixtures  were  found  to  be  inferior  to  the  CO^/N^/He  mixture  Single 
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pulse  peak  gain  of  the  CO^/N^/H.^  mixture  was  at  52%,  and  peak  gain  of  the 
CC^/CO/^/l^  combustion  mixture  was  at  32%  of  the  corresponding  value 
for  CC^/^/He  mixture  under  otherwise  identical  conditions. 

The  effects  on  performance  of  repetitively  pulsing  a recirculating  laser 
gas  mixture,  i.e.,  closed  cycle  operation,  were  investigated  for  the  CO^/CO/ 
N^/H^  combustion -generated  mixture.  It  was  shown  that  the  combination  of 
combustion  and  discharge -produced  contaminants  and  closed  cycle  operation 
may  severely  degrade  laser  performance.  It  was  also  shown  that  selective 
removal  of  the  most  important  electrophilic  reaction  products,  notably  HNO2 
and  HNOj,  promises  to  eliminate  the  problem  of  performance  degradation. 
Selective  removal  of  HNOx  can  readily  be  accomplished  between  passes  by 
passage  of  the  gases  over  a suitable  packed  bed  of  catalyst  (e.g.,  promoted 
dolomite)  since  just  downstream  of  the  cavity  the  gas  is  at  the  high  tempera- 
ture required  to  efficiently  use  this  process.  This  procedure  is  commonly 
employed  in  industrial  pollution  control  using  readily  available  technology 
and  commercial  catalysts. 

Because  of  the  constraints  imposed  on  the  development  of  the  present 
computer  code,  namely  maximum  use  of  existing  kinetic  data  and  of  a pre- 
viously developed  CW  electric  discharge  laser  code,  the  program  available 
now  is  highly  specialist -oriented  requiring  an  extremely  high  degree  of 
familiarity  on  the  side  of  the  user.  Applying  the  code  to  any  large  scale  or 
more  comprehensive  investigation  of  closed  cycle,  repetitively  pulsed  electric 
discharge  lasers  would  probably  be  both  impractical  and  uneconomical  unless 
some  effort  is  spent  on  refining,  streamlining  and  consolidating  the  program. 
Such  an  effort  should  aim  at:  yv)  the  further  refinement  of  the  present  code 
to  perform  detailed  analysis  of  certain  cases  of  particular  interest,  and  (b) 
developing  a "fast"  and  more  efficient  program  version  permitting  more 
global  analysis  of  the  essential  features  of  say  up  to  1000  repetitive  pulses 
by  sacrificing  a certain  amount  of  detail  available  from  the  "standard"  version. 
The  above  recommended  effort  should  most  definitely  include  further  develop- 
ment of  the  kinetics  model.  Formidable  complexities  and  uncertainties  are 
involved  in  modeling  the  plasma  chemistry  of  CC^  EDL  gas  mixtures  in  depth 
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sufficient  that  the  formation  of  electrophilic  species  concentrations  on  the 
order  of  parts  per  million  can  be  predicted.  In  addition  to  limitations  inherent 
in  the  information  gained  from  the  available  literature,  value  judgments  are 
required  with  regard  to  which  species  and  reactions  are  essential  to  incorporate 
in  the  model.  Minimization  of  the  number  of  species  and  reaction  mechanisms 
is  necessary  to  hold  down  the  required  program  core  storage  and  execution 
time.  The  latter  requirement  is  especially  severe  in  the  analysis  of  repeti- 
tively pulsed  lasers  for  which  computations  on  the  order  of  1000  pulses  might 
be  required. 

Confidence  in  the  logic  of  the  assumptions  inherent  in  the  program,  the 
need  for  revisions  or  extensions  and  the  degree  of  feasible  simplifications 
derives  most  productively  and  directly  from  further  comparison  of  theoretical 
and  experimental  results.  A valuable  experiment,  for  example,  would  be  to 
test  the  hypothesis  stated  in  Section  3 that  HNO^  might  be  one  of  the  culprits 
in  degrading  closed  cycle  pulsed  EDL  performance.  Such  an  experiment, 
using  a mass  spectrometer,  is  currently  planned  (Ref.  33).  Under  similar 
time-averaged  conditions  the  electron  concentration  is  necessarily  lower  by 
orders  of  magnitude  in  continuous  wave  operation  than  it  is  in  pulsed  operation. 
The  same  likely  holds  for  the  HNOx  concentrations  which,  under  CW  operating 
conditions,  may  be  so  small  as  to  be  undetectable  by  a mass  spectrometer. 

This  situation  points  to  the  importance  of  conducting  experiments  under  pulsed 
operating  conditions. 

Continued  data  analysis  and  correlation  would  greatly  enhance  the  signif- 
icance of  performance  predictions  and  as  well  provide  an  appreciably  greater 
insight  into  optimal  sustained  operation  of  closed  cycle,  repetitively  pulsed 
CO^  high  energy  lasers. 
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